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Abstract 



We have performed a detailed study of the zenith angle dependence of the 
regeneration factor and distributions of events at SNO and SK for different 
solutions of the solar neutrino problem. In particular, we discuss oscillatory 
behaviour and the synchronization effect in the distribution for the LMA solu- 
tion, the parametric peak for the LOW solution, etc. . . Physical interpretation 
of the effects is given. We suggest a new binning of events which emphasizes 
distinctive features of zenith angle distributions for the different solutions. We 
also find the correlations between the integrated day-night asymmetry and 
the rates of events in different zenith angle bins. Study of these correlations 
strengthens the identification power of the analysis. 
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I. INTRODUCTION 



Oscillations of solar neutrinos in the matter of the Earth modify the neutrino signal 
detected during the night [1-3] . The integral characteristic of this effect is the the day-night 
asymmetry: 

N-D 

, (1) 

where N and D are the night and the day event rates, averaged over the year, and corrected 
for the Earth orbit eccentricity. 

In most of the cases, the Earth-matter effect leads to the regeneration of the v e - flux, so 
the asymmetry is positive. Conversely, a negative asymmetry appears for the SMA solution 
at small mixing angles. There is a number of detailed studies of the asymmetry and its 
dependence on neutrino energy and oscillation parameters [1]- [14]. Asymmetries have also 
been studied for certain intervals of the zenith angles [6]. 

The observation of a day-night asymmetry will be the proof of the matter conversion 
solution of the solar neutrino problem, excluding the vacuum oscillation solutions. However, 
the measurement of the asymmetry alone may not allow to select among the three possible 
MSW solutions: LMA, SMA and LOW. A negative asymmetry is the evidence of the SMA 
solution. However the expected value of the effect is small |A DN | < 0.01 and it will be 
very difficult to establish experimentally such a small deviation from zero. Furthermore, 
one can make certain conclusions confronting the value of the asymmetry with other solar 
neutrino data (rates, energy spectrum of the recoil electrons). In particular, a relatively 
large asymmetry will be inconsistent with other solar neutrino data (rates, recoil electron 
energy spectrum) in the case of SMA and LOW solutions, thus favoring the LMA solution. 

Further insight can be obtained by studying the zenith angle distribution of events during 
the night. It turns out that different solutions lead to qualitatively different distributions. 
One of the first detailed studies of the zenith angle dependence of events for LMA, SMA 
and LOW solutions of the solar neutrino problem was done by Baltz and Wenesser [7] and 
then further elaborated in [4,8-11,14]. 

It was realized that for the SMA solution the distribution has a peculiar form with rather 
small effect for neutrinos crossing the mantle only (cos^z < 0.837) and large regeneration 
effect (peak) for neutrinos whose trajectories crossed both the mantle and the core of the 
Earth [7,9,10,15]. This peak at large cos^ was interpreted as due to certain constructive 
interference of the oscillation effects in the mantle and in the core of the Earth [16] or in 
more simple and transparent way, as the effect of parametric enhancement of oscillations 
[17] (see also [18,19] and for later discussion [20,21]). 

For the LMA solution one expects averaging of oscillations due to integration over the 
energies of neutrino and detected charge lepton as well as over finite size of the zenith angle 
bins. This leads to a rather flat distribution with small variations of the average rate during 
the night [7,11-13]. A significant regeneration effect is expected already in the first night 
bin cos#z < 0.2. 

For the LOW solution one gets the highest rate in the second night bin [7,14], followed 
by a dip and then again an increase of the rate for large cos^z- This peak has a simple 
interpretation as the oscillation maximum which corresponds to the phase of oscillations 

= 7T. 
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Thus, a detailed study of the zenith angle distribution will allow to disentangle the 
solutions. Moreover, it may allow to determine the oscillation parameters. Measurements of 
the zenith angle distribution are also interesting because they open the possibility to study 
various matter oscillation effects, such as: 

• oscillations in matter with constant and slowly changing density, 

• adiabatic conversion, 

• effects of the resonance enhancement of oscillations, 

• parametric effects, 

• oscillation effects in thin layers (with small density width), 

• effects of small density jumps, etc.. 

Furthermore, in principle, (if the oscillation parameters are determined from some other 
experiments) the precise measurement of the zenith angle distribution will allow to check 
models of the Earth density profile. 

In this paper we continue to study in detail the zenith angle distributions for the different 
MSW solutions of the solar neutrino problem. We give physical interpretation of various 
features of the distributions. We clarify what can be learned besides identification of the 
solution of the solar neutrino problem by measuring the zenith angle distributions in the 
present and in future high statistics experiments. 

The paper is organized as follows: in Sec. II we present general expressions for the 
survival probability and the regeneration factor which describe the Earth-matter effect and 
we study the properties of the regeneration factor in the constant and slowly changing 
density approximation. In Sec. Ill we present the results of our numerical calculations of 
the regeneration factor and the rates for the Super-Kamiokande and SNO for a realistic 
Earth-matter profile. We give an interpretation of the results of these calculations using 
the results of the analytical studies. In Sec. IV we propose a new cos 6z binning which 
emphasizes distinctive features of the distributions for the different solutions. We also study 
the correlations of rates as well as ratios of rates with the value of the day-night asymmetry 
as a way to strengthen the identification power of analysis. We summarize and discuss our 
results in Sec V. 



For the range of oscillation parameters of interest (Am 2 > 10 8 eV 2 ), the probability 
P ee to detect the solar electron neutrino at a detector can be written as follows: 



II. GENERAL RELATIONS 



A. Regeneration factor and conversion inside the Sun 




(2) 



In Eq.(2) the Earth regeneration factor, /reg, is defined as [14]: 



2 



/reg = P2e - sin 2 9 , (3) 

where P 2e is the probability of the v 2 — > v e conversion inside the Earth. In the absence of 
the Earth-matter effect we have P 2e = sin 2 9, so that /reg = 0. 

In Eq. (2) Pi is the probability of the v e — > u x conversion inside the Sun which can be 
approximated by the well known expression 

Pi = l+(l-Pc)cos29s . (4) 



2 V2 

Here 9s is the matter mixing angle at the production point inside the Sun: 

cos 29 s = cos 29 rn (p s ) (5) 
where 6 m {p) is the mixing angle in matter determined by 

nn — 1 + 7] cos 29 . . 

COS2 ^ = (l-2,cos2^ + , 2 )V 2 - ® 

Here 

Z V2m N Am 2 ( Am 2 / E Wig cm- 3 \ 

is the ratio between the refraction length, l , and the neutrino oscillation length in vacuum, 

lu. 

2nm N A-kE 

lo = V2^pV: l = a^- (8) 

In Eqs. (7) and (8) p is the matter density, Y e is the number of electrons per nucleon and 
ttl j\j is the nucleon mass. P c is the jump probability which takes for an exponential density 
profile the following form [22,23]: 



3 -7sin 2 6» _ -7 



1 - e-T ' 

where 7 is the ratio of the density scale height l p and the neutrino oscillation length: 



(9) 



-,- 4n% -i05( Am2/E )( l A l=-E- (10) 

r = i? o /10.54 is the height scale in the exponential approximation to the solar density 
profile. Inserting (4) into (2) we get 

P ee = P D -(l- 2P C ) • cos 29 s • ; re g , (11) 

where 

P D = ^ + ^(1-2P C ) -cos20 5 -cos20 (12) 
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is the survival probability in the absence of the Earth-matter effect, i.e., during the day. 
According to Eq. (11), the regeneration effect appears multiplied by two factors: 

1) The adiabaticity factor, (1 — 2P C ), which describes the adiabaticity of the conversion 
inside the Sun. The factor is maximal for P c = 0, that is, for the case of pure adiabatic 
propagation and it decreases with the increase of the adiabaticity breaking. For P c = 1/2 
the factor and the regeneration effect are zero, and at P c < 1/2 the Earth-matter effect 
changes the sign. 

2) The resonance position factor, cos2#s, which determines how far the resonance layer 
is situated from the production point (in the density scale). When neutrinos are produced 
at the resonance, one has cos2# 5 = 0, while cos 29 s < if the resonance occurs at densities 
lower than the one at the production point. The parameter cos2^5 decreases with the reso- 
nance density, and it reaches the value cos 29$ ~ — 1 when the resonance layer is sufficiently 
far out from the production region (i.e., at much lower densities). 

From the previous discussion we conclude that the largest regeneration effect occurs when 
the neutrino propagation is adiabatic and the resonance happens far from the production 
region. 

From Eq. (11) one can obtain the daily average survival probability which for rj <^ 1 
(cos29s ~ —1) takes the form: 

P = 1 -{P D + P N ) = ^[1 + (1 - 2P 1 )(/ reg - cos 29)]. (13) 

The numbers of events during the day and during the night are proportional to averaged 
probabilities Pd and Pjv, therefore the day-night asymmetry (1) can be written as 



MDN 



~ Pp 2/ re g 

P = l/(l-2P 1 )-cos2^ + / reg ' 



(14) 



where both in Eq.(13) and Eq.(14), / r eg should be averaged over the neutrino trajectories 
during the night. 



B. Regeneration factor in constant density and adiabatic approximations 

To obtain the zenith angle distributions we have performed exact numerical calculations 
of the regeneration factor integrating the evolution equation in the Earth-matter with the 
Earth density profile given in the Preliminary Reference Earth Model (PREM) [25]. How- 
ever, a number of qualitative features of our results can be easily understood in the simplified 
two-layers approximation of the Earth profile. In this approximation the profile consists of 
the mantle and the core with slowly changing densities and a sharp change of the density 
between the layers. In what follows we will parametrize the zenith angle dependence in 
terms of |cos#z| which, for simplicity, we will write omitting the moduli. Equivalently, 
this would correspond to 9z being the nadir angle. We call the region of zenith angles 
cos9z = — 0.837, for which neutrinos cross the mantle only, the mantle region and the 
region cos9z = 0.837 — 1, for which neutrinos cross both the mantle and the core, the core 
region. 

Simple analytical results can be obtained in the constant density approximation, or, 
in general, in the adiabatic approximation. In the case of constant density (which would 
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correspond to neutrinos crossing the mantle at small enough cos 6z ) we obtain the following 
expression for the regeneration factor 



2 i ird 



/ reg = sin2# m sin(2# m - 29) • sin 2 ^— J , (15) 

where 9 m is the mixing angle in the Earth matter given in Eq. (6) with 77 evaluated for the 
Earth density, and d is the width of the layer or equivalently the distance travelled by the 
neutrino. Using the expression for the mixing angle (6) we get that: 

77 sin 2 29 . 2 f nd\ 

freg = jz 7r—7r , ^ - sin — (16) 



3g (1 - 277 COS 20 + 77 2 ) 

which can be rewritten in the following form (see also [7]): 

/ reg = - ■ sin 2 26 m ■ sin 2 ( ?*) . (17) 

The oscillation length in matter, l m , equals 

sin2g m 77 k 

m v sin 2^ "(1 -2t/ cos 29 + r] 2 ) 1 / 2 (1 - 2t/ cos 20 + t? 2 ) 1 ^ ' ^ 

where Iq and / v are defined in Eq.(8) and 77 should be taken for the Earth density. 

Written in the form of Eq.(17), /reg differs from the expression for the probability of 
usual flavor oscillations in matter by the factor 1 / 77 oc n e E/Am 2 . One can notice as well 
that the amplitude of oscillations of / r eg 

A f = ^sin 2 29 m (19) 

is symmetric with respect to the exchange 77 <-> 1 / 77 , as can be seen explicitly from the 
Eq. (16). Due to the additional factor 1 / 77 the asymptotics of the oscillation amplitude are 
Af — > both at i] — > and at 77 — > 00 which differ from the ones for flavor oscillations [1]. 

According to Eq. (17), in the case of one layer with constant density, the regeneration 
/reg is always (for any value Am 2 / E and 9) positive, i.e. the matter effect of Earth always 
enhances the survival probability P ee (see Appendix). 

In fact, in the mantle (and also in the core) the density changes rather significantly. 
According to PREM model it increases from p ~ 3.2 g/cm 3 near the surface of the Earth 
to p ~ 5.6 g/cm 3 at the border with the core. The density changes smoothly apart from 
several jumps at distances (0.05 — 0.12)R E from the surface. Such a density variation leads to 
deviation from the simple oscillation picture describe above. One expects certain interplay 
of the oscillations and the adiabatic evolution which results in the change of the oscillation 
probability from the constant density description. Furthermore the small jumps in the 
density may induce some irregularities in the behavior of the zenith angle distribution. 

The description of the effects is simple if the adiabaticity condition is fulfilled (as it 
happens for the LMA solution). In this case the average probability and the amplitude of 
oscillations are determined uniquely by the instantaneous value of the density. Therefore P 
and the depth of oscillations of / r eg will be determined by the density at the surface p s : 
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/reg = -^r • sin 2 26 m {p s ) • sin 2 U f P) ■ (20) 

ViPs) V J ImJ 

In some cases when adiabaticity is violated one can describe the results using constant 
density approximation with some effective value of the density for each trajectory: 

l m = l m (e Z ), 0m = 0m(0z), V = V(0z)- (21) 

If the oscillations are averaged due to integration over the energy or/and the zenith angle 

we get 

f - ^ 20 (22) 

/reg 2(l-277cos2# + r? 2 ) 1 } 

which coincides with expression used in [14] with rj taking an average value for the corre- 
sponding E and 9 Z range. 

In the case of the core crossing trajectories the description becomes more complicated 
[16,17,20,21,30]. 

The analytical results presented in this section allow us to get a straightforward inter- 
pretation of the results of the numerical computations. Moreover, they give the correct 
functional dependence of observables on the neutrino parameters. 



III. ZENITH ANGLE DEPENDENCE OF THE REGENERATION FACTORS 

AND RATES 

In this section we consider the zenith angle dependence of the survival probability in 
Eq. (2), the regeneration factor in Eq. (3), and the event rates at SK and SNO for oscillation 
parameters from the LMA, LOW, and SMA solutions to the solar neutrino problem. 

The reduced rate of the electron scattering events [ES] at Super-Kamiokande is defined 

as 

[ES]" Nes 



/ dE' e R(E e , E' e ) f dE v F(E v )a e (E 1/ , E' e )[P ee (E u , 9 Z ) + r{E u , E' e )(l - P ee (E u , 6 Z )), (23) 
AW JE th J 

where R(E e , E' e ) is the energy resolution factor, F{E V ) is the flux of the electron neutrinos 
without oscillations, a e is the v e — e elastic scattering cross-section and r(E u ,E' e ) = a^/a e . 
with <t m being the — e elastic scattering cross-section [26]. N ve is the normalization factor 
which equals the integral in (23) taken at P = 1. 

The reduced rate of the charged current events [CC] at SNO is obtained from 

AT 1 

[CC] = JjStm = j Eth dE'MEe, K) j dE u F(E u )a cc (E u , E' e )P ee (E u , 9 Z ), (24) 

where acc is the CC neutrino-deuteron cross-section [27], and N u d is the normalization 
factor which equals the integral in (24) taken at P = 1. 
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The results presented in this section have been obtained for oscillation parameters in 
the presently allowed regions of solutions of the solar neutrino problem. In Fig. 1 we show 
the results of the global fit to the data which include (i) the SuperKamiokande (SK) data 
after 1117 days of operation (total number of events and day and night energy spectra), (ii) 
Gallex, GNO and SAGE data, (iii) the Homestake data. Shown are the best fit points for 
each solution as well as 90% and 99% CL regions found from the local minima in each allowed 
region. We show the solutions only in the range Am 2 > 10~ 8 eV 2 for which Earth-matter 
effects on the boron neutrinos can be substantial. The fit includes the latest standard solar 
model fluxes, BP00 model [28]. For details of the statistical analysis applied to the different 
observables we refer to Ref. [29]. 

The results for the event rates presented in Figs. 3-10 have been normalized to the 
SK total rate. That is, for each set of the oscillation parameters, expressions in Eqs. (23) 
and (24) have been multiplied by a boron flux normalization factor f$K m order get value 
of [ES] as it is measured at SK. 



A. LMA 

In the LMA region, the adiabaticity condition is satisfied and we can safely put P c = 0. 
Consequently Eq. (11) simplifies to: 

P ee = P D - cos 26 s freg, (25) 

where 

P ° = \ ~ \ C ° S 29s C ° S 29 ' ^ 
Moreover, in the low Am 2 part of the LMA region we can take cos 26$ ~ —1, so that 

P ee = sin 2 # + / re g, (Am 2 <10- 4 eV 2 ). (27) 

Let us discuss first the behaviour of the regeneration factor. In Figs. 3a and 4a we show 
the zenith angle dependence of the fceg f° r different values of neutrino energies. In the 
interval cos #2- = — 0.837 the dependence has a quasi-periodical shape. The amplitude of 
oscillation slightly changes with cos 9z as a consequence of the adiabaticity violation (arising 
from the small jumps of the density). At cosOz > 0.837 the dependence is more irregular due 
to the Earth core effect. All these features, as well as the dependence of the regeneration 
factor on the neutrino parameters (energy, mixing angle, and Am 2 ) can be immediately 
understood from our analytical consideration as we discuss next. 

In the LMA region of oscillation parameters the Earth-matter effects can be considered 
as neutrino oscillations in matter with slowly changing density. Moreover, in this region 
rj > 3 (for the present best fit point (see Fig. 3) and E ~ 10 MeV we get rj — 12). This 
means that oscillations in the Earth with parameters from the LMA region proceed in the 
vacuum dominated regime, when the matter effect gives relatively small corrections. In 
particular, the oscillation length is close to the vacuum oscillation length: 



/ ~ / 

<<m ~ l v 



2 cos 29 

1 + 

V 



(28) 
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which holds to first order in I/77 and presents a rather weak dependence on the mixing angle. 
For the best fit point oscillation parameters and pY e ~ 2 g/cc we get 

K 103 km (loSv) ' < 29) 

So, the oscillation length is much shorter than the diameter of the Earth for all the relevant 
oscillation parameters. 

The expression for the mixing angle can be approximated as 



tan26L tan 29 



1 + 



COS 20T) 

and the oscillation phase acquired by neutrinos crossing the Earth equals 



(30) 



^_ 2ird(6 z ) ^ 2nDcos9 z 

lm\Qz) K 

where D — 1.3 x 10 4 km is the diameter of the Earth. That is, the phase increases linearly 
with cos 9z, and therefore, /reg turns out to be an almost periodical function of cos 9 Z - Also, 
the phase is inversely proportional to the energy, so that the number of periods increases as 
1/E. The period of oscillations in the cos 0z scale equals T '(cos 9z) ~ l v /D. 

In fig. 2 we show the dependence of /reg 011 the distance travelled by the neutrino inside 
the Earth for some trajectories and for best fit values of the oscillation parameters. For 
LMA the oscillation length is much smaller than the typical scale of the density change: 

both for the mantle and for the core: l m <C l p . Since the mixing angle is large this leads to 
good adiabaticity (especially far from the resonance). Small density jumps can be treated 
then as small perturbations. As a result, the average /reg, /reg, follows the density change. 
Both / re g and the amplitude of oscillations increase towards the center of the trajectory 
where the density is maximal. Since the profile is symmetric (with respect to the middle 
point of the trajectory), the dependence of /reg and Af on distance is also symmetric. At 
the detector /reg and Af are determined by the surface density (p ~ 3.2 g/cc). Therefore 
in the adiabatic approximation the amplitude of the oscillatory behavior of /reg should not 
depend on cos 9 Z . The variations of the amplitude Af with cos 9z seen in Fig. 2 are produced 
by the small density jumps which violate the adiabaticity. 

As follows from Eq. (19), the amplitude of oscillations takes the form 

sin 2 29 1 , n . Ep s 
A , „ _ M e z ) „ _* , (33) 

where r] s and p s are the value of parameter r\ and the density at the surface of the Earth and 
k na ~ 1 is the parameter which describes small effects of the adiabaticity violation. Thus, 
the amplitude increases with energy, as can be also seen from Figs. 3a and 4a, and it is 
inversely proportional to Am 2 . The amplitude is proportional to sin 2 29, however variations 
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of this parameter in the LMA region produces small changes specially in the near maximal 
mixing region. 

Combining the results of Eqs. (31) and (33) we get the zenith angle dependence of the 
regeneration factor in the mantle region 



sin 29, .„ . . 9 tcD cos 9 z 
/reg « k na (0z) sin 2 . (34) 



The zenith angle dependence becomes more complicated in the region cos9z > 0.837, 
when neutrino trajectories cross the core of the Earth [16,17,20,21,30] . This is related to 
the appearance of the parametric (enhancement and suppression) effects on the top of the 
resonance enhancement in the core. Notice that for core-crossing trajectories the regenera- 
tion factor can be negative. In contrast, in the mantle region it is always positive as can be 
seen from Eq. (34) (see also Appendix). 

Finally let us point out that the results shown in Figs. 3- 7 for the regeneration factor 
are independent of the experimental set up. However the range of zenith angles covered 
by a given detector depends on its latitude: for SK cos9™ ax = 0.975, whereas for the SNO 
cos9™ ax = 0.92. For the sake of clarity in Figs. 3- 7 we mark those limiting zenith angles. 

Let us now consider the zenith angle dependence of the event rates. We have calculated 
the rates above a given energy threshold using Eqs. (23) and (24). The calculation of rates 
involves folding of the survival probability with neutrino cross-section, the flux of neutrinos 
and integration over the energy above threshold. Therefore results are different for SK and 
SNO. 

In the lower part of the LMA region the survival probability is simply given by the sum 
of sin 2 6* and the regeneration factor as seen in Eq. (27). Consequently, the [CC] rate can be 
written as 

[CC]^sin 2 ^+(/ re g(^))cc, (35) 

where the averaged regeneration factor is given by the equation (24) with P ee substituted 
by /reg- 

Similarly, the [ES] rate can be expressed as 

[ES]^sm 2 9 + rco S 2 9 + (l-r)(f Teg (9z)}ES , (36) 

where (/reg(^z)) es is given by Eq. (23) with P ee substituted by /reg- Notice that the 
interactions of and v T enhance the average [ES] rate in comparison with [CC] rate giving 
an additional term r cos 2 9, and they suppress the regeneration term by a factor (1 — r). 

In Figs. 3b and 4b we plot the rate [ES] at SK as a function of cos9 Z - The dependence 
of the reduced [CC] rate at SNO on cos 9 Z for two different thresholds is shown in Figs. 3c 
and 4c. (Similar figures have been obtained in [7,4]). As seen in the figures, the rates are 
oscillating functions of cos#z, however in contrast with the case of the regeneration factors, 
the amplitude of these oscillations changes with the zenith angle significantly. The amplitude 
is maximal at small cos9z {cos9z = — 0.2), it decreases with cos^z till cos 0z ~ 0.5, 
and then it increases again. The first maximum is achieved already at cos9 z = 0.02 — 
0.03 which corresponds to the distances 250 - 400 km. This large matter effect on small 
distances is related to the fact that the initial state is not a flavor state but an incoherent 



9 



admixture of the mass eigenstates v\ and [24]. Also the period of oscillations changes: it 
is about 0.07 for small cos9z and it decreases down to 0.05 for cos#z > 0.6. This important 
feature is related to the integration over the neutrino energy and it can be explained as the 
synchronization effect of the "waves" (oscillatory curves) /reg(^z) corresponding to different 
neutrino energies as we discuss next. 

Let us notice that the boron neutrino flux is maximal at E — 8 MeV and it decreases 
with energy, while both the cross-section and the regeneration factor increase with E. As a 
result, the main contribution to the integrated regeneration factor, (/reg(^z)), comes from 
a rather narrow interval of energies: E = 9 — 13 MeV, so that AE/E ~ 1/3. For the central 
value of the energy E = 11 MeV and oscillation parameters in the best fit point (see Fig. 3) 
the oscillation length is about 10 3 km. This distance corresponds to cos 9z ~ 0.07. That 
is, for neutrinos propagating at cos 9z ~ 0.07 the oscillation phase will be <fi ~ 2n and freg 
will be at a minimum. In general, we get from Eq. (31) 

, n „ „ ( A ™ 2 \ /10MeV\ 

= 147rcos0 z = — 9 — ^ — • 37 

r v 3xl ° eV A E ' 

Therefore, the phase difference A0 for neutrinos which cross the Earth in the direction cos#z 
and differ in energy by A E is 

. AE n cos9 z AE 

** K +E-"*'-0W-E- (38) 

If A0 ~ 2n one would expect strong averaging effect. From this condition and using Eq. (38) 
with AE/E = 1/3 we find that strong averaging should appear at cos9 z 0.2 in agreement 
with the results of the numerical calculations, as seen in Figs. 3b and 3c. For large cos9 z , 
the waves are again synchronized, so that a constructive interference of the waves leads to 
the restoration of the oscillatory behavior. 

From Eq.(38) we also see that the periods of constructive and destructive interference 
become shorter with the increase of Am 2 (see Fig. 4). Also with the increase of the energy 
threshold the relevant energy interval AE becomes narrower and, correspondingly, the period 
of synchronization becomes longer (in cos 0z scale) as can be seen in Fig. 3c and 4c. 



B. SMA 

In the SMA region we have cos 29 s ~ — 1- Furthermore, the adiabaticity is broken inside 
the Sun: P c ^ 0, so that 

P ee = P D + (l-2P c )-/ reg . (39) 

Notice that since in the SMA region cos 29 ~ 1, one has Pd ~ P c - SK data imply that 
P c ~ Pjj = 0.35-0.65, therefore in the SMA region the regeneration effect is substantially 
suppressed by the adiabaticity factor: |1 — 2P C | < 0.3. For P c > 1/2 the regeneration effect 
becomes negative. 

In Fig. 5a we show the zenith angle dependence of the regeneration factor for different 
values of neutrino energy and for oscillation parameters corresponding to the best fit point 
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in the SMA region. The behaviour of the curves can be understood taking into account that 
in the SMA region 



rj = 0.3 - 3 , (40) 

so the SMA solution region covers the resonance range for the Earth matter densities. For 
the best fit point we have 

l v « l « (4 - 8) x 10 3 km . (41) 

For the mantle region (cos9 z < 0.837) we find that: 

1. Far from the resonance, rj > 1, the oscillation length, 

* < 42 > 

is smaller than the Earth diameter. However the amplitude of oscillations is strongly sup- 
pressed (see the line which corresponds to the E = 5 MeV in fig. 5a) 

2. In the resonance, rj — 1, the amplitude of oscillations is strongly enhanced, but the 
oscillation length 

» — l T7TZ ~ 105 km > ( 43 ) 
sin 29 

is much larger that the Earth diameter. As a result, the oscillation effect is small (the line 
for E — 12 MeV in fig. 5). In the small phase limit the regeneration factor equals 

. ~ „ n 2 D 2 cos 2 9 Z , . , s 

/ reg » sin 2 29 z - . (44) 

It increases quadratically with o,o$9 z . 

3. The line for E — 15 MeV represents an intermediate case when the phase of oscillations 
is about 7r for cos#z ~ 0.8. In this case the regeneration factor equals: 



... D 2 cos 2 9 z , . L , D 2 cos 2 ^ 
/reg w 4 sin 2 29 — xs 4 sin 2 29 , (45) 



where in the second equality we have taken into account that r) ~ 1. Since Iq ~ 0.4D for 
large cos#z trajectories in the mantle, the mixing enhancement factor, 4(_D// ) 2 , can be as 
large as 25, so that / r eg ~ 0.04 for sin 2 29 = 2.4 x 10~ 3 and cos^ z = 0.8. 

Clearly, the adiabaticity is strongly broken near the resonance: l m ^> l p . But far from the 
resonance: l m ~ l p ~ D, the violation of the adiabaticity is moderate, so one can describe 
its effect as oscillations in a narrow layer with some effective density. 

For cos 9z > 0.837 neutrinos cross the core of the Earth. For the core densities the 
resonance energies are in the range E = 3 — 5 MeV (Am 2 ~ 5 x 10 -6 eV 2 ). For energies 
between E R (core) and E R (mantle), the parametric enhancement of oscillations takes place 
leading to the appearance of the parametric peak in the cos#z distribution of /reg- For 
SMA the realization of the parametric resonance corresponds to a mixing angle in mantle 
smaller than maximal mixing, 29 mant i e < 7r/2, while in the core, 29 core > n/2. At the 



11 



peak the relation between the oscillation phases in the mantle 4> ma ntie and in the core cor e 
corresponds to the general condition for the parametric resonance [17,20,21] 

A 3 = s m c c cos 29 mantle + s c c m cos 29 core « , (46) 

where s m = sin(0 mant / e /2), c c = cos(0 man ^ e /2). For the best fit point tan 2 9 = 6.1 x 10 -4 , 
Am 2 = 5 x 10 -6 eV 2 and E = 10 MeV we find that the maximum of the parametric peak 
occurs at cos9 z = 0.863. The phases in the core and in the mantle are 4> ma ntie = 0.49-7T, 
(ficore = 0.497T, and the effective mixing angles are 29 man tie = 0.188 and n — 29 core = 0.113. 
Then we get X 3 ~ 0.045 < 1 . 

As illustration, in Fig. 2b we show the dependence of the regeneration factor on the 
distance for the parameters which correspond to the maxima of the parametric peaks. 

In Figs. 5a and 6a we show the dependence of the regeneration factor on the mixing 
angle. Both the oscillation length and the resonance condition depend on sin 2 29 weakly, so 
that with a good precision / r eg oc sin 2 29. Notice that this proportionality holds also for 
core crossing trajectories. 

In Figs. 5 and 6 (panels b and c) we show the zenith angle dependence of the event rates 
at SK and SNO. The integration over energy basically reproduces the dependence of /reg on 
cos9 z for E ~ 10 MeV since this energy gives the dominant contribution. The [CC] event 
rate at SNO can be writen as 

[CC] » (P D ) CC + (1 - 2P c )(f reg (9 z ))cc , (47) 

where P c is an effective jump probability for the contributing energies. (Notice that P c 
changes slower with E than /reg does.) 

For tan 2 9 ~ 6 x 10" 4 (see Fig. 5) the adiabaticity breaking is strong, P c > 1/2 and the 
negative adiabatic factor leads to the suppression of the rate due to Earth-matter effect. 
On the other hand, for tan 2 9 ~ 2.0 x 10~ 3 (see fig. 6), the adiabaticity breaking is weaker, 
P c < 1/2, and the Earth-matter effect is positive. As follows from the figures, the adiabatic 
factor suppresses substantially the regeneration effect: 

|l-2P c | (/reg f z))cC <0.2-0.3. (48) 
/reg 

The [CC] rate depends weakly on the energy threshold in the mantle region but it strongly 
decreases with E t n in the core region. The parametric peak is wider (in energy scale) than 
the resonance peak in the mantle and therefore the integration over the energy leads to a 
stronger decrease of the regeneration effect in the mantle region than in the core region. 
Also the parametric peak is situated at the low energy part of the spectrum: (5 - 10 MeV). 
As a consequence, the core peak decreases substantially with the increase of the threshold 
from 5 to 8 MeV. 

C. LOW 

In the LOW solution region one has cos2#s ~ — 1, and moreover, to good approximation 
we can take P c ~ 0, so that 
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P ee « sin 2 9 + /reg- 



(49) 



The zenith angle dependence of the regeneration factor for oscillation parameters in 
the best fit point and different values of the neutrino energy are shown in Fig. 7a. The 
interpretation of the results is rather straightforward. In the LOW region we have 77 < 0.1. In 
particular, in the best fit point at E = 10 MeV : rj < 0.03. Therefore the oscillations proceed 
in the matter dominated regime. Thus, l v ^> Iq and the oscillation length is determined 
mainly by the refraction length. Notice that in the limit C 1 we get from Eq. (18) 



l m « Z (l + 2r/cos29) 



(50) 



Moreover, in the LOW region the mixing parameter is small: cos 29 < 0.5 at 99% C.L. (in 
the best fit point: cos 29 ~ 0.2). So, the correction to Iq is further suppressed, for instance, 
2?7cos2# ~ 0.012 for E — 10 MeV. Therefore for a given trajectory, the oscillation phase 
is practically independent of the neutrino energy and mixing angle (see different curves in 
Fig. 7a). 

The phase of the oscillatory behavior of /reg with cos 9 Z can be written as 



rDcos8 z 

4>=y2Gp I n e (x,cos9z)dx 
Jo 



2irD cos 9 Z 
l (cos9 z ) 



(51) 



where l (cos9z) is the average refraction length along the trajectory determined by cos 9z 
and n e = pY e /m N is the electron number density. With increase of cos 9z, the average 
density increases so, l decreases and the period of oscillations becomes shorter. Thud the 
oscillation phase is determined by cos9z- The first maximum of /reg (0 = 7r ) is achieved at 
cos 9z = 0.35, the minimum (0 = 2n) lies at cos 9z = 0.59 and second maximum occurs at 
cos#z = 0.8 (see also [7]). 

The oscillation length is comparable with the scale of density variations l m ~ Iq ~ 8 • 10 3 
km, where l p ~ D/2 , so that the adiabaticity is broken moderately. Since l m is comparable 
with the size of the layer, the effect can be considered as oscillation in the layer of matter 
with some effective constant density. 

The amplitude of oscillations: 



A f ~ rj sin 2 29 



v^Am 2 sin 29 
G F En(cos9 z ) 



(52) 



is proportional to Am 2 and inversely proportional to E. It is also inversely proportional to 
the average density for a given trajectory n which increases with cos^z. Correspondingly 
the second peak is lower. For E= 10 MeV and oscillation parameters at the best fit point 
we get rj = 0.2 (for the surface density), and consequently according to Eq. (52) the height 
of the first peak is expected to be /reg = 0.042 which is slightly larger than the numerical 
result /reg = 0.035 shown in Fig. 7a. The difference of the results is due to the adiabaticity 
breaking. 

Combining Eqs. (51) and (52) we get an approximate expression for the regeneration 
factor in the mantle region: 



/] 



V2 sin 2 29 Am 2 



reg 



Gpn(cos9z) E 



sin 



Q F rD cos 9 Z 
V2J0 



n e (x, cos 9z)dx 



(53) 
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Notice also that /reg increases with sin 2 29. 

In the core region / re g is enhanced due to the parametric effect. Indeed, we find that 
for cos 9z = 0.92, which corresponds to the position of the third maximum, the phases of 
oscillations equal: 

<Pmantle = 0.98vT, (f) core = 2.98tT . (54) 

That is, they are very close to 7r and 3n 

^mantle ~ 7T, 0core ~ 37T . (55) 

At this condition in the constant density approximation the height of the peak would reach 
the value (see Appendix): 

f Te g ~ sin 2 (6> - A9 mantle + 29 core ) - sin 2 9 m sin 29 ■ sm 2 (A9 mant i e - 26» core ), (56) 

where 9 ma ntie and 9 core are the mixing angles in the mantle and in the core correspondingly, 
and the last equality is valid for 9 man tie ~ 9 core ~ n/2. The peak in Fig. 7a is slightly 
lower than what Eq. (56) gives because of the density change. Notice that this realization 
of the parametric enhancement corresponds to both mixing angles above the resonance, 
(see [30] and Appendix). The real time evolution of the neutrino state (dependence of the 
regeneration factor on distance) for the parameters corresponding to the peak is shown in 
fig. 2c. 

In Fig. 7b and Fig. 7c we show the zenith angle dependence of [ES] event rates and [CC] 
event rates at SK and SNO respectively. The [CC] rate has the same approximate expression 
as in Eq. (35) for the LMA solution. Since / re g ocl/E and the phase of oscillations depends 
very weakly on E, we get (freg)cc ~ freg(E), where E ~ (10 — 11) MeV is the effective 
energy of the spectrum. Therefore the [CC] distribution reproduces all the features of the 
zenith angle dependence of /reg- Furthermore, the regeneration term weakly decreases with 
the increase of the threshold energy from 5 to 8 MeV, and it increases with mixing as 
{freg)cc sin 2 29 according to Eq. (53). 

For the [ES] distribution at SK, the results are similar to those for [CC] with the only 
additional feature of a 25 % damping effect due to the contribution from v^j v T scattering 
via NC (see approximate formula in (36)). 

IV. IDENTIFICATION OF THE SOLUTION. CORRELATIONS 

A. Binned rates 

According to Figs. 3- 7 the LMA, SMA and LOW solutions present qualitatively different 
zenith angle dependence of the event rates (either [ES] or [CC]). 

1) For the LMA solution one expects a significant rate already at cos9 z ~ 0.03. The first 
peak is at cos9 z = 0.02 — 0.03. The rate has an oscillatory behavior at small, cos#z < 0.2, 
and large cos#z > 0.6 — 0.8, zenith angles. The border of the oscillatory region at high 9z 
depends on the value of the mass difference. For instance, for Am 2 ~ 2.5 x 10~ 5 eV 2 the 
second oscillatory region starts at cos9 z ~ 0.8 while for Am 2 ~ 4.5 x 10~ 5 eV 2 it expands 
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down to cosOz ~ 0.6. The period of oscillations is small, A(cos#z) < 0.05 — 0.07, and in 
consequence, it will be difficult to detect the oscillatory behaviour due to the relatively low 
statistics. Although we find that for small values of Am 2 the amplitude of oscillations can 
reach 15 % of the average rate. When averaging over wider cos9 z bins, the binned rate 
depends rather weakly on cos 9 Z and no significant change of the rate is expected in the core 
region. 

2) For the SMA solution the rate changes slowly and monotonously in the mantle re- 
gion and significant parametric effects can appear in the core region. For the large tan 2 9 
the rate increases with cos9 z and the parametric enhancement leads to the appearance of 
the parametric peak which gives the main contribution to the integral regeneration effect 
[7,8,10,4,5]. For small tan 2 9 the rate decreases with cos9 z and in the core one has a dip in 
the distribution of events. 

3) The LOW solution predicts the existence of three peaks in the zenith angle distribution 
of events. There are two oscillation peaks in the mantle range: a wide peak with the 
maximum at cos9 z = 0.35 and a narrow peak at cos9 z ~ 0.8. The third peak (in the core 
region) is due to the parametric enhancement of oscillations. The maximum of this peak is 
at cos9 z ~ 0.92. 

These features are rather generic. The qualitative behaviour of the distribution is the 
same for all points within a given solution, although quantitatively the sizes of the different 
structures appearing in the distributions change with the oscillation parameters. Therefore, 
in principle a detailed study of the zenith angle distributions will not only give the identifi- 
cation of the solution of the problem but also the determination of the neutrino oscillation 
parameters. 

Due to the relatively small present statistics (effects themselves are small) it is unavoid- 
able to bin the distribution of events. The binning should be chosen in such a way to avoid 
as much as possible the averaging or washing out of the structures. 

We propose the following binning. We will enumerate the night bins as Nl ... N5 and 
denote by [N]i, (i = 1...5) the average rate in a given bin Ni normalized to the no oscillation 
prediction. We will call it the reduced rate. We will denote by [N] and [D] the averaged 
reduced rates during the night and during the day correspondingly. 

Since in the mantle range both the SMA and LMA solutions predict a rather flat dis- 
tribution with weak dependence on cos9 z , we suggest the binning which emphasizes the 
features of the distribution for the LOW solution. 

• Bin Nl: cos9 z = 0-0.173 (90° > 9 Z > 80°). In this bin one expects very small 
regeneration effect both for the SMA and LOW solutions. Significant effect should be 
observed in the case of the LMA solution. 

• Bin N2: cos9 z = 0.173-0.5 (80° > 9 Z > 60°). This bin is selected in such a way that 
it covers the main part of the wide peak of the LOW distribution. 

• Bin N3: cos0 z = 0.5-0.707 (60° > 9 Z > 45°). This bin covers the dip of the LOW 
distribution: [N] 3 (LOW) ~ [D]. 

• Bin N4: cos9 z = 0.707-0.83 (45° > 9 Z > 35°) is the last mantle bin. This bin 
corresponds to the second peak of the LOW solution. 
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• Bin N5: cosdz > 0.83 (9z < 35°), the core bin. This bin is restricted by cosOz = 
0.83-0.975 for SK and cos9 z = 0.83-0.92 for SNO. For LMA the rate is comparable 
with the rates in the previous bins. For the LOW solution one expects to see a rate 
slightly smaller than the one in the fourth bin. Here the regeneration effect due to 
the parametric enhancement of oscillations takes place. For SMA the effect depends 
strongly on the oscillation parameters: for large tan 2 9 one should see the highest rate. 
With decrease of the mixing angle the rate decreases and becomes smaller that the 
day rate. 

In Figs. 8, 9, and 10 we show the binned event rates for [ES] events at SK and [CC] events 
at SNO for oscillation parameters in the LMA, SMA and LOW solutions respectively. The 
distributions are rather similar for [ES] events at SK and [CC] events at SNO. Although the 
enhancement of the Earth-matter effect in the core bin for the SMA solution is larger in 
SK, since this bin covers a larger interval of cos 9z at the SK latitude. 

The qualitative behaviour of the [CC] rates does not change with the increase of the 
threshold energy from 5 MeV to 8 MeV for SNO. The absolute value of the effect (reduced 
rate), however, increases for the LMA solution and it decreases for the LOW and SMA 
solutions. 

For the LMA solution the regeneration effect decreases with the increase of Am 2 in all 
the bins and it changes weakly with the mixing. For the LOW solution the situation is 
opposite: the rates [N]2, [N]4 and [N]s increase with Am 2 . This behaviour originates from 
the dependence of / r eg on Am 2 (see Eqs. (53) and (56)). Also for SMA the rate in the core 
bin increases with tan 2 9. 

The comparison of the event rates rates in the different bins allows to discriminate among 
the solutions of the solar neutrino problem. In particular, for LMA we have 

[N] x < [N] 2 < [N] 3 < [N] 4 < [N] 5 . (57) 

For LOW distribution: 

[N] 2 > [N] 4 > [N] x ~ [N] 3 > [D]. (58) 

For SMA distribution: 

IfN], - [D]| < |[N] 2 - [D]| < |[N] 3 - [D]| < |[N] 4 - [D]| < |[N] 5 - [D]|. (59) 

Finally, for the sake of completeness we show in Fig. 11 the time exposure for the different 
bins (average during the year) for the SK and SNO detectors. The histograms are normalized 
in such a way that the integral over night time gives 0.5. The time exposure distribution 
determines the statistics in each bin. As follows from the figure the highest statistics is 
expected in N2 followed by N3. Moreover, at SNO the time exposure N2 will be about twice 
longer than N3, while for SK the difference in statistics is smaller. For SNO the lowest 
statistics is expected in N5, whereas for SK it corresponds to N4. 

B. Correlations 

At present, the expected regeneration effects are restricted by the Super-Kamiokande 
result on the day-night asymmetry: 
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(60) 



The asymmetry gives the integrated (over cos6z) Earth-matter effect. With present statis- 
tics the integral effect is at 1.3<r level which means that in spite of large number of events 
it will be difficult to measure the zenith angle distribution with a high precision. Higher 
asymmetry (due to absence of damping) is expected for SNO 



A^P < 0.10-0.15. (61) 

To enhance the identification power of the analysis one can study the correlations between 
the signals in different zenith angle bins and the integral effect which can be represented by 
the day-night asymmetry. As an illustration of such correlations we first show in Fig. 12 con- 
tours of constant asymmetry Adn and of ratios of the reduced rates [N] 2 /[N] 3 and [N] 5 /[N] 2 
in the Am 2 - tan 2 9 plane for SNO with threshold 5 MeV. 

For the LMA solution the cos 6z - distribution is rather flat which means that for all bins 



m_ 1 



oc A DN . (62) 



[N] 

For the LOW solution the heights of all three peaks increase with integral effect: 

[N]s 



[N] 



- 1 oc A DN (63) 



and similar proportionality exists for [N]4 and [N]s. In the third bin, however, the regener- 
ation effect is practically zero and the rate is at the level of the day rate. 

For the SMA solution, the rate monotonously increases (or decreases depending on tan 9) 
with cos9 z , so that for all the bins the proportionality (62) holds approximately. 



1. Ratio [N] 5 /[N] 2 versus Aqn- 
Comparing the signals in N2 and N5 bins we get 

[N]5 - ■ 1 (LMA), & > 1 (LOW) . (64) 



[N] 2 v " [N] s 

In both these cases the day-night asymmetry is positive. On the other hand for SMA we 
have [N] 5 /[N] 2 > 1 for the part of the allowed region with small mixing, for < 0, and 
[N] 5 /[N] 2 < 1 for the part of the region with large mixing, where A^n > (see Fig. 12). 

In Figs. 13, and 14 we show the results of mapping the 99% CL solution regions onto the 
plane of the [N]s/[N] 2 — Adn observables for the SNO and SK experiments. From the figures 
we see that there is a strong correlation between the [N] 5 /[N] 2 ratio and the asymmetry for 
all three solutions. 

For SNO observables above 5 MeV (Fig. 13) the correlation can be parametrized as: 



[N] ; 



1 LMA 

1 - 0.45A DN LOW . (65) 
1 + 3.6A DN SMA 
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A negative asymmetry would testify for the SMA solution but the maximal asymmetry 
which is expected for the SMA on the basis of the present data cannot be larger than 0.05. 
Thus, for A DN (SNO) > 0.05 (above 5 MeV) one will choose between LMA and LOW. A 
ratio [N] 5 /[N] 2 > 1 will testify for LMA, whereas [N] 5 /[N] 2 < 1 will be an evidence of the 
LOW solution. For A DN = 0, the equality [N] 5 /[N] 2 = 1 is fulfilled for all the solutions. In 
the range A D n(SNO) < 0.05 all three solutions are possible. With increase of the asymmetry 
the difference of ratios for different solutions increases. 

For the SNO measurements with E t h = 8 MeV, the correlations are described approxi- 
mately by (65), but the LMA region has a more irregular shape since at E th = 8 MeV the 
averaging is weaker. The maximal allowed deviations of [N] 5 /[N] 2 from 1 are ~ 11% (SMA), 
~ 4% (LMA) and ~ -5% (LOW). 

Comparing with the results at lower energy threshold we find that the allowed SMA 
region in [N] 5 /[N] 2 — A DN plane shrinks and A DN (SNO) can reach only 0.04, while [N] 5 /[N] 2 
can still be as large as 1.11. This happens because the region of parametric enhancement 
is less covered with the higher threshold. For the LOW solution the allowed region in- 
creases (since the regeneration effect decreases with energy). Adn(SNO) can reach 11% and 
(N 5 /N 2 - 1) ~ 0.05. 

For the ES - events at SK (Fig. 14) one gets a similar picture, although the regeneration 
effects are damped by the and u T contribution. The expected maximal asymmetry is 
A DN (SK) ~ 3% for SMA and A DN (SK) ~ 10% for LMA and LOW. If A DN (SK) ~ 5%, the 
[N] 5 /[N] 2 ratio for LMA is larger than the ratio for LOW by (0.03 - 0.05). 

To have an idea about present sensitivity in searches of correlations, we also show in 
Fig. 14 the SK results and with 1 a errors for the day-night asymmetry and the [N]s/[N] 2 
ratio. One must bare in mind, however, that the plotted result on [N] 5 /[N] 2 does not 
correspond to the binning we are proposing. In Fig. 13 we show the corresponding attainable 
SNO sensitivity for three years of operation (corresponding to 13200 CC events for 5 MeV 
threshold and 6000 CC events for 8 MeV threshold). For definiteness we have plotted the 
point at a value Adn(SNO) = 1.6Adn(SK) which is the expected relation for the best fit 
point in LMA as discussed in Ref. [32]. 



[N] 2 » [N] 3 (LMA), [N] 2 > [N] 3 (LOW), |[N] 2 - [D]| < |[N] 3 - [D]| (SMA) . (66) 



So that the rates in the high statistics bins N2 and N3 can provide an important criteria to 
distinguish the solutions. 

In Fig. 15 we show the mapping the 99% CL allowed regions of oscillation parameters 
onto the [N] 2 /[N] 3 — A D n plane of the SNO observables for two different energy thresholds. 
From the figure one can see that there is a clear correlation between the ratio and the 
asymmetry for the SMA and LOW solutions: 



2. Ratio [N] 2 /[N] 3 versus A DN . 



According to Eqs. (57)-(59) (see also Fig. 12), we have 




LMA 
LOW 
SMA 



(67) 
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Also, we see that the regions only weakly depend on the energy threshold. In the case of 
the LOW solution the [N] 2 /[N] 3 - ratio can reach 1.17 for E th = 5 MeV and 1.15 for E th = 8 
MeV. For SMA the maximal allowed deviation of the [N] 2 /[N] 3 ratio from 1 is about 4%. We 
get qualitatively similar results for SK (Fig. 16), where the regeneration effect is damped by 

and v T contributions, so that the allowed regions for all three solutions are smaller than 
for SNO. We also plot in Figs. 16 and 15 the present results from SK and the estimated 
accuracy at SNO as discussed at the end of the previous section. 

From these results it is clear that, with the estimated statistics at SK and SNO, it will 
be difficult to get strong discrimination among the solutions. At 3<r level all the solutions 
seem to be allowed. However, future higher statistics experiments like UNO [33] can give a 
more significant result. 

One may also attempt to elaborate a test criteria for some of the solutions without a 
large lost of statistics. For instance, for the LOW solution one expects a significant deviation 
from one for the ratio: 

[N] 2 + [N] 4 + [N] 5 

[NK + [N] 3 > (68) 

While for LMA this ratio should be close to one. 



V. CONCLUSIONS 

In this work we have studied in detail the zenith angle dependence of the regeneration 
factor and the distributions of the solar neutrino events expected in the Super-Kamiokande 
and SNO detectors. 

We have identified the generic features of the distributions for each solution. 

1. The LMA cos6z - distribution is characterized by a fast oscillatory behavior with 
rather slowly changing averaged regeneration rate. The amplitude of oscillations is enhanced 
in the core region. The behaviour of the event rates at SNO and SK has an interesting 
"synchronization effect" so that strong oscillations of the rates with the change of cos $z are 
only expected at small and at large values of cos 6z- In particular 2-3 oscillations with 
significant amplitude are expected at small cos 9z- 

To detect such a behavior one needs to use small size binning A(cos9z) ~ 0.05 which (in 
view of relatively low statistics) may only be possible with new megaton water Cherenkov 
detectors. If the LMA solution is identified, a study of the oscillatory behavior can be used 
to measure the oscillation parameters (in particular Am 2 ) and the Earth density profile. 

2. The SMA solution predicts a smooth zenith angle dependence in the mantle region 
with small regeneration effect while the parametric effects can show up in the core region. 
For the larger mixing part of the allowed region the regeneration effect is positive with a 
parametric peak in the core region; for the smaller mixings the effect is negative and in the 
core region one should see the parametric dip. 

3. In the case of the LOW solution there are two oscillation peaks in the mantle region 
with maxima at cos 9z = 0.35 (oscillation phase <f> = ir) and at cos 9z = 0.78 (the phase 3n). 
In the core region the third peak appears due to parametric enhancement of oscillations. In 
the maximum the conditions for the parametric resonance are approximately satisfied. This 
realization of the parametric resonance corresponds to the case of neutrino energy above the 
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resonance energy in the mantle (which differs from the realization for the SMA) and to the 
phase relations <fi ma ntie ~ tt and core ~ 3-7T. We find that the position of the peaks is almos 
independent of the neutrino energy, whereas the height of the peaks as well as the integral 
effect are inversely proportional to the energy. 

The shape of the distributions and their dependence on the oscillation parameters as well 
as energy thresholds are well described by oscillations in a medium which consists of one 
or three layers with slowly changing density. Such an oscillation picture gives the correct 
functional dependence of the rates on the neutrino parameters. 

We also show that precise measurements of the zenith angle distribution will allow 
to study a number of various matter effects such as oscillation in matter with constant 
and slowly changing density, adiabatic conversion, parametric enhancement of oscillations, 
etc 

We have suggested a new binning in the zenith angle distribution which is chosen to 
emphasize the distinctive features of the distributions for the different solutions. The iden- 
tification power of the analysis can be further enhanced by studying correlations between 
the event rates in the different bins and the integral effect which can be represented by 
the day-night asymmetry. In particular we have shown that the study of the correlation 
between the ratio [N] 2 /[N] 3 and the day-night asymmetry will help to identify the LOW 
solution while the correlation between [N] 5 /[N] 2 and A D / N is more suitable for identification 
of the SMA solution. The distributions expected at SNO and at SK have similar properties. 
But at SNO the structures in the distributions and the asymmetry are enhanced due to the 
absence of the damping effect. 
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APPENDIX. GRAPHIC REPRESENTATION OF EVOLUTION OF THE 

NEUTRINO STATE 

We describe here a graphic representation of the evolution of the neutrino state which 
allows us to understand the properties of the zenith angle distributions for the different 
solutions. The graphic representation [19] is based on the analogy of the neutrino evolution 
with the behaviour of spin of the electron in the magnetic field. The neutrino state is 
described by a vector of length 1/2 with components 



where ipi, (i = e, /i) are the neutrino wave functions of the electron and muon neutrinos. 
(The components of this vector are elements of the density matrix.) 
Introducing the vector 
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(69) 



B = — (sin20 m , 0, cos2# m ) 



(70) 
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(6 m is the mixing angle and l m is the oscillation length in the medium) which would corre- 
spond to the magnetic field, one can get from the evolution equation for ipi the equation for 
v: 

£=(*xa). (7i) 

According to (69) the projection of V on the axis z, u z , gives the probability to find u e 
in the state V: 

P = ^ = v s + \ =cos 2 |. (72) 

Here v z = O.Scos^e, and 6 e is the angle between z7 and the axis z. The z axis can be called 
the flavor axis as it coincides with vectors z7 e = — z7 M . 

In a medium with constant density, 9 m = const, the evolution corresponds to the preces- 
sion of the vector V - around B: v moves according to the increase of the oscillation phase, 
0, on the surface of the cone with axis B. The vector B coincides with the vector of the 
eigenstates in matter z/ im = — z7 2rrt . We will denote by V2 m [mantle) and z7 2m (core) the axis in 
the mantle and in the core. Notice that when the matter density is below the resonant one 
(77 > 1) the corresponding axis lies in the first-third quadrants with decreasing projection 
on the z axis as the density increases (see Figs. 17a and 17b). Once the resonant density is 
crossed, the axis moves to the second-fourth quadrants (see Figs. 17c and 17d). 

In Figs. 17a and Figs. 17b we show the evolution of the neutrino state in the case of the 
LMA solution for trajectories crossing the mantle only (a) or the mantle and the core (b). 
Inside the Sun the neutrino vector (produced as V = u e ps z7 2m follows the vector z7 2m due to 
adiabaticity and appears at the surface of the Sun and then at the surface of the Earth as 
v = i?2- Due to the adiabatic density change in the mantle the axis i? 2m moves accordingly, 
but at the detector it (as well as the cone of precession) will return to the same position 
as in the moment of arrival of the neutrino at the Earth. The cones with axis z7 2m and P' 2m 
describe oscillations at the surface of the Earth and in the central part of the trajectory 
correspondingly. 

At the surface of the Earth the density jumps suddenly and therefore the axis vector 
U2m suddenly changes its position to from z7 2 to U2m(niantle) . In the mantle the neutrino 
vector will precess around z? 2m (mantle). Clearly, since V2m{jnantle) is closer than z7 2 to v> e , 
the rotation will lead to an increase of the z projection and therefore to increase the P ee 
probability. This illustrates the fact that the regeneration is always positive in the mantle 
zone (see also [31]). The quasi-adiabatic movements of z? 2m (man£Ze) due to density variations 
in the mantle do not change this conclusion. 

Let us now consider the evolution for the core crossing trajectory. The neutrino vector 
starts at z7 a = u 2 and it will first precess around u 2m (mantle). The phase of oscillations is 
rather large ~ (10 — 20)ir and it depends sensitively on the neutrino energy and the zenith 
angle. So, depending on the energy the neutrino state can arrive at the core in any position 
on the cone. In Fig. 17b, as an example, we have selected the state v^. In the core the vector 
V will precess around ^ 2m (core) starting from the position z7&. The phase is large and the 
state will enter the mantle again in some position V c which depends on E and Q z . In the 
second mantle layer the state will precess around z7 2m [mantle] starting from V c state. Let us 
denote by Vf the final neutrino state at the detector. As it is clear from the figure, in most 
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of configurations determined by the positions of V c and uj, (and therefore E and 6 Z ) the final 
state will have the projection u z larger that the projection of u a = v 2 (which corresponds 
to day signal). In other words: cos 29 e (iJf) < cos26>. In some cases, however, one can get 
v z ( final) < u 2z which corresponds to negative regeneration. 

It is easy to see that maximal regeneration effect would correspond to zero regeneration 
effect in the mantle: (fimantie = 2nk and to maximal effect in the core: core = n + Ink' (k, k' 
are integer numbers). 

In Fig. 17c we show the evolution of the neutrino state in the case of SMA solution. 
Now the neutrino arrives at the Earth as an incoherent superposition of u x and v 2 . This 
split of the state originates from the adiabaticity violation inside the Sun. The total survival 
probability P ee is then determined by independent oscillations of v\ and v 2 inside the Earth. 
The result can be represented in terms on v 2 — > v e probability and the adiabaticity factor 
(12). The later can have both negative and positive signs depending on the oscillation 
parameters, and it is this factor which determines the sign of the regeneration effect. 

In Fig. 17c we show the evolution of the v 2 state inside the Earth. In the mantle the 
vector V precesses around u 2rn {mantle) starting from V a = v 2 . In the core region, the 
parametric enhancement of oscillations occurs for certain values of cos#z- We show the 
graphic representation of evolution for this case. Neutrino arrives at the core in the state 
Ub, then V precesses around V 2m (core) starting from V h and reaches the border of core and 
mantle in the state V c (which corresponds to the phase cor e ~ 0.757r). In the second mantle 
layer V rotates around P 2m (mantle) starting from v c and arrive at the detector in the state 
Vf. Clearly, the regeneration turns out to be enhanced. 

In the case of the LOW solution (Fig. 17d), v arrives at the Earth as V a pa v 2 similarly 
to the LMA case but in contrast to LMA, now the mixing is matter suppressed so that 
the axis of the eigenstates is close to the flavor axis (z). In the mantle z7 precesses around 
V2m,{jnantle) starting from position V a . The peaks in the zenith angle distribution in the 
mantle region correspond to the final position of the neutrino v^. In the core region the 
parametric enhancement of oscillations occurs. The maximum of the parametric peak cor- 
responds to the following picture. Vector V precesses in the mantle and reaches the state 
V h at the border of the core (4> ma ntie ~ vr). In the core V rotates around u 2m (core) starting 
from z/fe. The neutrino vector make 1.5 turns (4> cor e — 37r) and leave the core in the position 
v c . In the second mantle layer V rotates again around V 2m {mantle) with initial position z7 c . 
At the detector one will detect the state which corresponds to Pf. 
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FIGURES 




FIG. 1. Regions of solutions of the solar neutrino problem in the Am 2 — tan 2 9 plane. Dots 
correspond to the best fit points for each type of the solution. Contours show 90 % CL regions 
(light) and 99 % CL regions (dark) found with respect to local minima for each region. 
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FIG. 2. The dependence of the regeneration factor on distance in units of the total length of a 
trajectory for a given value of cos 9z and for the oscillation parameters indicated in the pannels. 
The dependence is shown for (a) LMA solution, (b) SMA solution for two different energies (arrow 
and triangles indicate the point when neutrino enters and leave the core), (c) LOW solution. 
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FIG. 3. The zenith angle distributions for the best fit point (Am 2 = 3.25 x 10~ 5 eV 2 , 
tan 2 9 = 0.365) in the LMA solution, (a) The dependence of the regeneration factor on cos 9z 
for different neutrino energies as labeled in the figure. Vertical dashed lines indicate the maximal 
values of cosOz which can be realized at SK and SNO. (b) The zenith angle dependence of the ve 
event rate at Super-Kamiokande above the energy E t ^ = 5.5 MeV. (c) The zenith angle depen- 
dence of the vd charged current event rate at SNO for two different energy thresholds (5 and 8 
MeV ). 
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FIG. 4. Same as Fig. 3 for a larger Am 2 ( Am 2 = 4.5 x 1CT 5 eV 2 ). 
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FIG. 5. Same as Fig. 3 for the best fit point ( Am 2 = 5.01 x 1CT 6 eV 2 , tan 2 6 = 6.1 x 1CT 4 ) 
in the SMA solution. 
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FIG. 6. Same as Fig. 5 for a larger mixing (tan 2 9 = 2.0 x 10" 3 ) in the SMA solution 
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FIG. 7. Same as Fig. 3 for the best fit ( Am 2 = 1.0 x 1(T 7 eV 2 , tan 2 9 = 0.668 ) in the LOW 
solution. 
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FIG. 8. The binned zenith angle dependence of (a) [ES] event rates at SK and (b) the vd CC - 
event rate at SNO above at Eth = 5 MeV and (c) Eth = 8 MeV in the LMA region. The histograms 
correspond to different values of the oscillation parameters as labeled in the figure. The first bin 
corresponds to the day rate. The event rates have been computed with a boron flux normalization 
chosen to fit the measured rate at SK. 
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FIG. 9. Same as Fig. 8 for the SMA solution. 
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FIG. 10. Same as Fig. 8 for the LOW solution. 
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FIG. 11. The exposure time for the different cos9z bins for SNO and SK detectors during the 
year. 
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FIG. 12. The contours of equal Day-Night asymmetry, Adn (a), the ratios of rates N5/N2 (b), 
and N 2 /N 3 (c) at SNO ( threshold 5 MeV ) in the Am 2 - tan 2 9 plane. 
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FIG. 13. The allowed values for the day-night asymmetry, j4dn and the ratio of rates in the 
N5 and N2 bins, N 5 /N 2 , at SNO above 5 MeV (a) and 8 MeV (b) in the different allowed regions 
(at 99% CL) of the solar neutrino problem. We also show the expected statistical sensitivity after 
3 years of operation (see text for details). 
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FIG. 14. The allowed values for the day-night asymmetry, j4dn and the ratio of rates in the 
N5 and N2 bins, N5/N2, at SK in the different allowed regions (at 99% CL) of the solar neutrino 
problem. The present experimental values and 1 a errors are also displayed (see text for details). 
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FIG. 15. The allowed values for the day-night asymmetry, Adn and the ratio of rates in the 
N2 and N3 bins, N 2 /N 3 , at SNO above 5 MeV (a) and 8 MeV (b) in the different allowed regions 
(at 99% CL) of the solar neutrino problem. We also show the expected statistical sensitivity after 
3 years of operation (see text for details). 
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FIG. 16. The allowed values for the day-night asymmetry, j4dn and the ratio of rates in the 
N3 and N2 bins, N3/N2, at SK in the different allowed regions (at 99% CL) of the solar neutrino 
problem. The present experimental values and 1 a errors are also displayed (see text for details). 
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FIG. 17. Graphic representation of the evolution of the neutrino state inside the Earth for 
different solutions of the solar neutrino problem, (a) LMA solution: propagation along trajectory 
which crosses the mantle only; (b) LMA solution: evolution for core crossing trajectory; (c) SMA 
solution: evolution for core crossing trajectory and cos 9z which corresponds to the maximum of 
the parametric peak; (d) the same as in (c) for LOW solution. 
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